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The analysis of electromagnetic propagation over uniform cables 
depends on the calculation of the charge densities on the conductors, 
relative to a potential function that is not necessarily constant on the 
conductors. By considering such a potential function as the real part 
of an analytic function, two Laurent series are derived, one of which 
involves the Fourier components of the potential function and its as- 
sociated charge densities on the conductors. The second series accounts 
for the relative location of the conductors. The two series are equated 
to give a system of linear equations that can be solved for the charge 
components. The results obtained, which apply to uniform cables whose 
conductors (including the shield, if present) have circular cross sections 
and are covered with two layers of dielectric insulation, can be used to 
calculate the propagation modes and propagation constants of the 
cable. 

I. INTRODUCTION 

Multiconductor cables for telecommunications have distributions of 
charge on each conductor. The surface charge density on the conductors 
is proportional to the normal derivative of a potential function (i.e., a 
solution to Laplace's equation), which is defined in the region separating 
the conductors and is constant on the conductors. The constant values 
of the potential are the voltages of the conductors, and the proportion- 
ality constant is the permittivity of the material next to the conductor. 
Also, generalized charge densities are defined for potential functions, 
such as the longitudinal component of the electric field, which are not 
constant on the conductors, 

In the present work an algorithm is developed for computing the 
charge densities in this generalized sense for uniform cables. Thus, the 
conductors of the cable are assumed to be straight and parallel, so that 
each transverse cross section is identical. The wires are assumed to have 
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Fig. 1 — Typical cable cross section 

circular cross sections and to be covered by two circularly symmetric 
layers of homogeneous dielectric material. Surrounding the collection 
of wires is a circular, metallic shield, and it too is assumed to have two 
uniform layers of dielectric on its inside surface. A typical cross section 
is shown in Fig. 1. 

Interest in the charge densities was spurred by the recent finding that 
the modes of a cable and their associated propagation constants can be 
expressed directly in terms of the charge densities when, as in Kuznet- 
sov's work, 1 low frequencies are excluded. Subsequently, it was found 
that the analysis could be extended to low frequencies by using the notion 
of generalized charge densities associated with the longitudinal com- 
ponent of the electric field. These matters are developed in detail in 
Ref. 2. 

Previous work on calculating charge on the conductors of a cable 
usually involved the simplifying assumption of homogeneous dielectric 
in the region separating the conductors. For two identical wires in free 
space, an explicit formula is available for all the Fourier components of 
the charge density on a wire (see Ref. 1, page 41). Goluzin, 3 using the 
theory of complex variables, has developed a continued fraction ex- 
pansion for a potential function in a region bounded by circles; the ex- 
pansion converges under certain conditions on the size and location of 
the circles. Nordgard 4 has developed an algorithm involving a matrix 
inversion for computing Fourier components of the charge densities for 
a pair in a shield. Also, the capacitance elements (the zero-order com- 
ponent of a charge density) have been calculated in a variety of cir- 
cumstances. 5,6 

For inhomogeneous dielectrics, the calculation of charge density has 
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been approached by a finite-difference technique 7 coupled with a Fourier 
analysis 8 of the normal derivative to give the Fourier components. This 
technique, by calculating more information than is needed, has proved 
to involve more computer expense than the technique to be described 
here. 

In this paper, techniques from the theory of complex variables are used 
to develop a set of linear equations for the charge components. Following 
Goluzin (see Ref. 3), the potential function is viewed as the real part of 
an analytic function; then two Laurent series are calculated for it about 
each conductor. The first is expressed in terms of the Fourier compo- 
nents of the charge density and the potential function on each conductor; 
the second is expressed in terms of singularities located at the centers 
of the conductors (at infinity for the shield). Since the two Laurent series 
must be identical, their corresponding coefficients can be equated to give 
a system of linear equations from which the charge components can be 
calculated. For convenience, the system of equations is expressed in 
matrix form. 

The matrix equations are developed in the next section, with details 
of the calculations and proofs relegated to appendices. In Section III, 
the development is summarized and then extended to apply to cables 
without shields, to cables with holes in the dielectric, and to more general 
boundary problems that arise in determining the admittance matrix and 
other propagation parameters of the cable. Numerical examples and 
experimental testing of this work is presented in Ref. 9. 

II. CALCULATION OF THE CHARGE COMPONENTS 

The cable to be considered consists of M straight and parallel wires 
with circular cross sections surrounded by a circular shield of inside ra- 
dius Rn. The wires and the inside surface of the shield are covered by two 
different layers of dielectric insulation. For the mth wire (1 ^ m s M) 
the radius is R m ; the permittivity of its first layer of dielectric is t m with 
thickness R m \ ~ R m ; and the permittivity of its second layer is t m \ with 
thickness R m 2 ~ R m i- On the shield (referred to as the m — conductor), 
the permittivities of the first and second layers are c and eoi with 
thickness R - R ol and R m - Rm, respectively. The permittivity of the 
material separating the insulated conductors is e. 

When the cross-sectional plane of the cable is viewed as the complex 
plane, the centers of the conductors can be specified as complex numbers. 
These are denoted b m (m = 0, 1, • • ■ , M), where 6 = refers to the center 
of the shield. A typical cross section is shown in Fig. 1. 

With {p m , m ) as polar coordinates based at the center of the mth 
conductor (m = 0, 1, ■ • • , M), a potential function U is assumed to satisfy 
the following conditions: 

(i) In the region separating the conductors, V 2 C/ = 0, where V 2 = 
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(d 2 /dx 2 ) + (d 2 /dy 2 ) and (x, y) denote Cartesian coordinates in the plane 
(i.e., U satisfies Laplace's equation). 

(«) At the surface of the mth conductor (p m =R m ),U has the Fourier 
series: 

[/ = Eu«m«pN,n) m = 0, 1,---,M. (1) 

n 

(Hi) When e denotes the permittivity as a function of position, then 
both U and e(dt//dp m ) are continuous across dielectric interfaces about 
the mth conductor for m = 0, 1, - ■ • , M. The charge density associated 
with U at the mth conductor is 

dU 

Pm ~ € m ~ 

dPm 



= E Pnm exp (m<f>m) m = 0, 1, • - ■ , M. (2) 

Thus, the problem is to determine the Fourier components, p nm . 

The boundary problem for U is put into the context of complex vari- 
ables by viewing the cross-sectional plane of the cable as the complex 
plane. Then, in the region separating the conductors, there is an analytic 
function f(z) (unique within addition of an imaginary constant) such 
that U is the real part of f(z), And f{z) can be represented as a Laurent 
series. 10 

In Appendix A, a Laurent series is calculated for/U) in a neighborhood 
of b m in terms of the Fourier series of U and € m (dU/dp m ) at the surface 
of the mth conductor [eq. (1) and eq. (2), respectively] for m - 0, 1, 
* • • , M. The result, eq. (26) of Appendix A, shows the coefficients of the 
Laurent series about z = b m depending only on p nm and u nm (n = 0, ±1, 
±2, ■ ■ •) with no explicit indication of interconductor coupling. 

A second representation for f{z), which is based on Cauchy's integral 
formula is 



f(z)=f (z)+ E fm(*h O) 

m=l 



where 



A)(Z) = 0OO+ £ Pno(^-Y (4) 

is analytic everywhere inside the shield and 

f-nM-^en^^+id^^^y" (5) 

is analytic everywhere outside the mth wire for m = 1, • ■ ■ , M (see Ref. 
3). The coefficients j3 o, ■ • • , 0OM are real, and in general the /3 nm (for 1 
£ m ^ M and n > 0) are complex. In Appendix B, a second Laurent se- 
ries is calculated for /(z) in a neighborhood of b m by combining these 
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functions. The result is eq. (38) for m = and eq. (39) for m = 1, • • • , 
M. 

The two Laurent series must be identical since they represent the same 
function f{z). Therefore, their coefficients can be equated, and, as shown 
in Appendix C, this leads to systems of linear equations in u nm , p nm , and 
nm for m = 0, 1, - - • , M and n = 0, 1, - - • . It suffices to deal with non- 
negative values of n because U and p m are real quantities for all m 
(hence, u- nm = u* nm and p- nm ■ p* nm , where * denotes complex con- 
jugate). 

The systems of equations are conveniently expressed in matrix form 
with the various coefficients combined into vectors. Accordingly, the 
following infinite vectors are defined: 



(6) 





U m = (W(> m , U\ m , • ■ •) 




Pm = \Pl)m, P\m, ' ' *) 




ftn = WQm> 01m, • ' •), 


all for m - 0, 1 , - ■ 


• , M and the joint vectors, 




u = (uo, * - - , u m ) 




P = (PO, • • • , Pm) 




= (ft). ■■■ ,^m)- 



(7) 

As indicated in Appendix C, when these vectors are taken as column 
vectors, there are matrices T, G, and H such that 

T0 = u and p = Gp - Hu. (8) 

In particular, G and H are such that 

p 0m = (t/R m )0 Om m = l f - . . , M. (9) 

If only the zero-order components of the charge-densities are of in- 
terest, then it suffices to invert the T-matrix to give 

0=T~ l u, (10) 

whereupon eq. (9) is used. When all the components of the charge den- 
sities are of interest, then combining the equations in (8) gives 

p = (GT~ 1 -H)u. (11) 

(That the T-matrix is invertible follows from the well-known uniqueness 
of solutions to Laplace's equation with prescribed values on the 
boundary.) 

In practice, the infinite vectors u m , p m , and fi m m = 0, 1, • • • , M are 
truncated to give Af-vectors, and the matrices T, G, and H are truncated 
to (M + l)N X (M + 1)N matrices. The matrix operations indicated in 
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eq. (11) are then carried out on the truncated matrices to give an ap- 
proximation to the first N components of the charge density on each 
conductor. 

A detailed study of the matrix truncation has not been carried out, 
though the cables studied in Ref. 9 provide some experience in this 
matter. For the cables whose wires and shield were mutually separated 
by more than one wire diameter (the 754E and Focal), it more than 
sufficed to consider five harmonics on the conductors, including the zero 
order (i.e., N = 5). When N was set to 6, there was a difference in the 
coefficient of the dominant harmonic (zero order) of less than 1 in 10,000; 
furthermore, the coefficient of the extra harmonic was four orders of 
magnitude less than that of the zero-order harmonic. For the cable whose 
wires and shield were separated by a small fraction of one wire diameter 
(the Proximity cable), eight harmonics were required: the coefficient 
of the seventh-order harmonic was about % 5 that of the zero-order har- 
monic. When N was set to 10, again there was a difference in the domi- 
nant coefficients of less than 1 in 10,000; and the extra coefficients were 
two orders of magnitude less than that of the zero-order harmonic. 

A second practical matter is the presence of conjugation operators in 
the T-matrix. When it is known (e.g. , by symmetry) that the coefficients 
$ im are real for all £ and m, then the conjugation operator has no effect 
and can be ignored. The case where the conductor centers are collinear 
is handled in this way. 

In general, the coefficients have an imaginary part. Then $ (m is ex- 
pressed as the two-vector 

/Rej3 fm \ 

\Im/W 

and the element of the T-matrix, T km {n t £), which multiplies p em , is 
expressed as 

/Re Tkm (n, t) -* m T km <«. t)\ 

\JmT km (n,t) ReT km (n,£)/ 
The first component of the matrix product is the proper real part and 
the second is the proper imaginary part of the product, Tk m ( n > £)fiem- 
When a conjugation operator appears, the 2 X 2 matrix above must be 
multiplied on the right by 



C -5 



(this corresponds to changing the sign of the imaginary part of (m ). The 
result is the matrix 

/Re T km (n, £) Im T km (n, t) \ 
\lmT km {n,£) -BjeT km (n,t))' 
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When these matrices are substituted for the elements of the T-matrix, 
then the resulting matrix can be inverted in the usual way. 

III. SUMMARY AND EXTENSIONS 

Matrix equations have been derived which relate u, p, and the auxi- 
liary vector £. These are 

Tp = u and p = Gp - Hu. 
From these equations, the vector p can be determined when u is given. 
When the potential function U is the constant 1 on the mth conductor 
and zero on the others, the quantity 

C km = 2wR k pQk = 2tt€/3 * 1 ^ K m ^ M (12) 

is the fern- element of the M X M capacitance matrix C. Thus, the ca- 
pacitance matrix is calculated as a special case of the analysis in Section 
II. 

These results can be readily extended to cables without a shield and 
with some modification to cables with circular holes in the interstitial 
dielectric. When no shield is present, two changes must be made: 

(i) Components of u, p, and (3 associated with m = must be elimi- 
nated, and the corresponding submatrices in the matrices T, G, and H 
must be eliminated. 

(ii) The T-matrix must be bordered by one row and one column to 
give 

f= ( c , £) withc = ( e i'-"-' e i>> < 13 > 

where ei = (1, 0, 0, - - -) is repeated M-times. Then the first equation in 
(8) becomes 

'©-£)■ 

where Q is the total charge of the cable and £ is a constant to be deter- 
mined. The second condition comes from the requirement that the total 
charge be specified, but that the potential on the boundary (the surface 
of the conductors) be specified only within an additive constant (£). This 
holds for any exterior problem for Laplace's equation in two dimen- 
sions. 11 

Circular holes in the interstitial dielectric are treated like extra con- 
ductors with unspecified potential values at their surface (p/, = R/ } ). But 
if u n t t denote the Fourier components of the potential function at the 
surface of the hole, then inside the hole, 

— ) exp(in0/,); (15) 

tifj/ 

UNIFORM CABLES— II 617 



and if e h is the permittivity of the material in the hole, then 

dU 



out OPh 



= L {nthUnh/Rh) exp (in0,,), (16) 



dU 

Oph 

where the normal-derivatives are evaluated on the outside and inside 
surfaces of the hole as indicated. Therefore, in terms of their Fourier 
components 

V h =D h u h , (17) 

where Dh is the diagonal matrix with main diagonal \nthiRh } for n = 0, 
1, 2, • • - . Also, by eq. (46) 

p h = G k h - H h u h . (18) 

From eq. (17) and eq. (18), h can be calculated in terms of u h . This 
in conjunction with the matrix-equation T0 = u is sufficient to determine 
first fa, then 0, and then by the second equation in (8), p. Details are not 
supplied here. 

The problem can also be generalized by specifying more complicated 
boundary conditions on the conductors. For example, the longitudinal 
component of the electric field satisfies a boundary condition of the 

form, 

p - Svl = g, (19) 

(see Ref. 2) where S is some matrix and g some vector. It follows imme- 
diately from (8) that 

(G - HT -ST)0 = g. (20) 

Thus, under certain conditions on matrix S, is obtained by inverting 
(G -HT - ST) t u is obtained as T0, and p is then obtained from eq. 

(19). 

The latter problem is involved in determining the elements of the 
admittance matrix for the cable (see Ref. 2, eq. 21). As shown in Ref. 2, 
this is an intermediate step for calculating the propagation modes of the 
cable and their associated propagation constants. The finite-difference 
technique, which had been used to calculate charge densities, 7 - 8 was not 
capable of dealing with this type of boundary condition. But even for 
cases such as the calculation of the capacitance matrix where the finite 
difference technique could be applied, the technique described here has 
a cost savings of more than one order of magnitude for 0.1 percent ac- 
curacy. Thus, the technique has proved to be flexible in solving potential 
problems associated with uniform cables and relatively inexpensive. 

APPENDIX A 

7fce First Laurent Series 

When the cross-sectional plane of the cable is viewed as the complex 
plane, the potential function U is the real part of an analytic function 
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f{z). In this appendix a Laurent series is deduced for f(z) in a neigh- 
borhood of each conductor, starting from the Fourier series for the charge 
density and U at the surface of the conductor, as given in eq. (1) and eq. 

(21- 

In the first layer of insulation for the mth conductor (for R m = p m = 
fffli m = 1, ■ • • , M or i? = Po = Rqi), the potential function is 

U(Pm, <Pm) = U 0m + p 0m {R m /t m )£n £m + % £ 

X exp (m<p m )[u nm ($% + £-") + Pnm(RmMt m )(F m - £-")], (21) 

where £ m = (p m /R m ). This is validated by noting that it satisfies La- 
place's equation and it satisfies the boundary conditions of eq, (1) and 
eq. (2) at p m = R m (i.e., when £ m = 1). 

In the second layer of insulation (for fl ml s p m ^ R m2 m = 1, • • • , M 
or R oi ^ po g R02)* the potential function is 

U(pm, 0m) = "Om + P0m(Rm/t m ) ^ nr m + Pom (R m /*m l)^"Ul 



+ y 4 E exp(in0 m )[u nm G;L(p m )+p nm (fl m /n« m )G nm (p m )], (22) 
where £ ml - (p m /R m \), 

r m ■ (Rml/Rm), 5m = (€ml/€m) (23) 

and 

G£Jp m )-(ri±r;")tt5 1 , + E;!) 

+ (*J- 1 (r:,Tr-")(r:,-U). (24) 

This satisfies Laplace's equation and the continuity conditions on £/ and 
HdU/dp m ) at the interface p m ■ R m L (i.e., when f m i — 1 and £ m ■ r m ), 
For t/ this is obvious, but to check the continuity of l{bUf&p m ), it is 
useful to refer to the calculation, 

inli> m )H$ m ip m )m — G$ m ip m ) = (n/p m )[(r n m ±r~ n )(f ml - £"?) 

Spm 

+ («J- 1 (rlTr„7)(r;, + U)]. (25) 

Outside the insulation (for p m ^R m2 +,m = 1, - - * , M or p < R02), 
the analytic function whose real part matches U is 



f{z) = Uom + POmiRm/thm + Pi\ m {Rmh)tn ( ~ ~ ) 

\ R m / 

00 /z — b \ n 
+ V4 L ( '") [wnmA nm +p nm (/Z m /ne m )B nm ] 

+ y 4 £ f^i- 21 ) "(«;^ fl m + p; m (fl m /n€ m )F Bin j, (26) 
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where * denotes complex conjugation, 

K m ■ &mi&mtnr m + 6 m i£nr m] - £nr m r m i, (27) 

r m \ - (Rm2/Rml), 5ml = («/<ml), t 28 ) 

and with G* m = G* m (i? m2 ) and //„„, = H* m (R m 2), 

A nm = (r m r ml )-HGt m + i^mi^Htm) 

Bam - (r m r ml )-"(G- m + ($mi)~ l ffnJ 

£ nm = (r m r m i)"(GL - (*ml) _1 HJJ 

F rtm - (r m r ml nG- m - (5 m i)- x tf,7j. (29) 

Substituting z = b m + R m2 exp (intf> m ) into f(z) and taking the real part 
gives 

U() m +Pom(^mA)(«m + fWm'mi) 

+ % E exp (in0 m )[u nm (-A„ m (r m r m] ) n + ^n m (r m r m] )~ rt ) 

+ p„ m (fi m /nc m )(B rtm (r m r ml )" + F nm (r m r m i)- n )], 

since u_ nm = u* m and p~ nm = pL- As is easily checked, this matches 
U in eq. (22) for Pm = R m2 (i.e., for Ui = r m ]l Likewise, it is straight- 
forward to check that the p m -derivative of the real part of f{z) matches 
Umi/* m )dU/dp m evaluated at Pm = R m2 -. This validates eq. (26) for 

Hz). 

When the dielectric is homogeneous, then r m - r mi - d m - b m i - 1. 
It follows that in this case K m = and A m = B nm = E nm = —F nm = 4. 



APPENDIX B 

The Second Laurent Series 

An analytic function f(z) in the region separating the conductors has 
the form given in eq. (3) through eq. (5). In this appendix, a single Lau- 
rent series is derived for f(z) in a neighborhood of each conductor (in- 
cluding the shield) by combining these equations. 

The functions / m (z), as given in eq. (4) and eq. (5), are analytic in a 
neighborhood of b k for m ^h and k 9* 0; hence, they are represented by 
a power series about z = bk, 

f m (z) = £ c m k(n) (^r^Y m = 0, 1, . • • , M. (30) 

n=0 \ tik f 
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The coefficients c mh (n) are related to the derivatives of f m {z) at z = b k 
by the formula 

c m k(n)=—*- — ^f „ =0,1,2,.--; (31) 

nl dz" \z=bk 

in particular, c mh (0) = f ni (b b ). 

The results of the calculation in eq. (31) are indicated at the end of 
this appendix. When the formula is evaluated for f n (z), as given in eq. 
(4), the result is eq. (40), and for f m (z), as given in eq. (5), the result is 
the combination of eq. (41) and eq. (42). The combined Laurent series 
for f(z) about z = b k is indicated in eq. (39). 

In a neighborhood of the surrounding shield, the Laurent series of eq. 
(5) holds for f m (z) (m = 1, • • - , M). This is rewritten as 



where 



fmiz) = £ 0m in (~j + (3 0m en ^ + g m {z), 



(32) 



(33) 



Since g m ( °° ) = 0, it follows that g m {Rl/z ) is analytic in a neighborhood 
of z - 0; so it is represented there by a power series, 

'-(f)".?, 8 -^©"' (34) 

and accordingly, 

gm(z) = L c m0 (n) (-£-) ". (35) 

The coefficients c m0 (n) for m = 1, - . • , M are obtained from the for- 
mula, 



, v Rl\ d" f /Rhl 



2=0 



n = 1,2,-... (36) 



The result of applying this formula to eq. (33) is 

n = 1, 2, - • • 

The combined Laurent series for f(z) in a neighborhood of the shield 
is 

f{z)= i(3no(f-Y+ E \dmin(~-) 

+ A»m ** (|^) + £; Cm0 („) (^)~"1 (38) 
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The combined Laurent series for f(z) about z = b k (h ^ 0) is 

iz — bb\ ~ n /z — bh\~ n 

f(z) = j£ w(0) + a» m (— ) + e a* (^-) 

The coefficients obtained by carrying out the calculations in eq. (31) 



are 



and form,i2^0,m^fe 



c mk 



c m k(0) = £ 0m £n 

/6 m -6ft\-« 



b fc - 6, 



R, 



fcy; - &m\ ' 



+ £ <w( « ) 



(40) 



(41) 



<"' = (^) 



"[-^m'*: -1 )^)"] <42) 



for n = 1, 2, 



APPENDIX C 

Equate the Two Laurent Series 

Equating the constant, logarithmic, nth power, and -nth power terms 
(n ^ 1), respectively, in the Laurent series about the shield [i.e., in eq. 
(26) and eq. (38)] gives the linear equations, 

■ # n 

(l) 000+ E 0mO*"-^- = "OO + Poo(flo/«)Ko 

m = l *im 

M 

(») E 0om = Poo(fioA) 

m = l 



<«l) p\i0 = y^UnO^nO + Pno(^o/"eo)Sno] 



M 



(iu) E c m o(n) - %["« Ao + pBo(^o/««o)^no]- 



(43) 



m = l 



For the fcth wire (k = 1, • • • , M), the corresponding equations are 
(i) E c mfe (0) = Ko* + Pofc(fi*A)«* 

(ii) /3 ft - poh(Rhfi) 

(Hi) E c mfe (n)= y 4 [u n A^„fc + p n A(R*/n«fc)BnA] 

(w) fink = Uu nk E nk + p nk (RkIntk)F nk l (44) 
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The charge coefficients can be solved from {ii) and {Si) in the first set 
and from {ii) and {iu) in the second set to give 

M 

POO = U/Ro) E &Qm 
m = l 

Pno — (4eo"/#oJ3no)#nO _ UonA n o/RoB n0 )u n0 

p Qk = U/Rk)f$ok 

p nk ■ (4t k n/RkF n k)P' n k ~ UknE nk /RkF n k)unk- (45) 

In terms of the vectors defined in (6) (viewed as column vectors), these 
equations are expressed 

P* = G k (3 k -H k u k k « 0, 1, • • • , M, (46) 

where the n^-elements of these matrices for n t £ — 0, 1, -•• are 



Gain, £) = 



(0; ei,***; ei) n = 
|(4e n/fioB n o) n = i * 
1 otherwise 



G k (n,£)= U/Rk) n = £ = Q 

(k= l.-.-.Af) |(4£ A n/i?^„ A )(*) n =/ ^0 
n ?£ ^ 



Ho(n, /) = 
(*-l,---,M) 



{f {) nA f]0 /^oB n o) n = £ 

n ^ / 

\U k nE nk /R k F nk ) n = £ 

n j* £ 



In terms of the joint vectors defined in (7), 

p = G0 - Hu, (47) 

where G and H are the direct sum of the matrices, G k and H kt respec- 
tively, for k = 0, 1, • • • , M (see Ref. 12, p. 159). 

The charge components can be eliminated in the original set of 
equations by means of the identity, 

BnkEnh — A n hF n k = (32/5* 5* 1) 
n = 0, 1, • • • and k = 0, ■ * ■ , M. (48) 
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The result is the set of equations, 

(0 0oo + E 0om ( in - KQ J = uoo 

M 
m=l 

("0 £ C m fc(0) - 0OfcKfe = UQk 

(U)) - (&khlFnk/B) £ C mk (n) + (dkhlBnk/W*nk = Unk 

for fc = 1, - - • , M and n - 1, 2, • • • . 

These equations can be expressed in the matrix form 



Tfi = 



TooToi • • • Tqm 
TiqT u - •• Tiji/ 



00 




ruoi 


01 


= 


Ui 


l^ m - 




_UM_ 



= u, 



(49) 



(50) 



TmqTmi * * * Tmjw_ 
where and u are viewed as column vectors. The symbols Tkm for £, m 
= 0, • • • , M denote submatrices of T defined as follows: 

n = £ = 



T 00 {n,£) = 



1° 



n ^ t 






Tomin, t) = 
(m = 1, • * • , M) 



= ^ = 






n = 1,2,--- 

otherwise 



(^'„» 0i / = o,i,... 

\ . /A /bk\~ n /°ft\' l^n<^ 



7Wn, /) = 



/ = 1, 2, • • 

otherwise 
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T k k(n, £) = 
(k = l,...,M) 



•Kk 



n = £ =0 



(hihBnk/BK*) n = £ * 

n^t 



b m ~ b k \ "" t = 0; 

Rk I n = 1,2, • • • 

i + n - 1\ /fr m - b fe \ -" /6 fe - 6,-\ 



T* m (n f /)= Inl 6 * H n=* = 

(fe, m^O Jfe^m) ( bfe ~ m ) ' n = 0;£ = l,2, 

\ R m / 

«, ^ ^ o. 

The symbol (*) throughout means that complex conjugation is to be 
performed. These elements of the various submatrices are read off the 
appropriate coefficients in the system of equations in (49) in conjunction 
with the equations for the c mk {-) in eqs. (37), (40), (41), and (42). 

REFERENCES 

1. P. I. Kuznetsov and R. L. Stratonovich, The Propagation of Electromagnetic Waves 

in Multiconductor Transmission Lines, New York: Pergamon Press, 1964, Chap 
1,2. 

2. T. A. Lenahan, "The Theory of Uniform Cables— Part I: Calculation of Propagation 

Parameters," B.S.T.J., this issue, pp. 597-610. 

3. G. M. Goluzin, "Solution of the Principal Plane Problems of Mathematical Physics 

for the Laplace Equation and Multiply Connected Domains Limited by Circles 
(Method of Functional Equations)," Mat Sbornik, 41, No. 2 (1934), pp. 246-276. 

4. J. D. Nordgard, "The Capacitances and Surface -Charge Distributions of a Shielded 

Balanced Pair," IEEE Trans. Microw. Theory Tech., 24, No. 2 (Feb. 1976), pp 
94-100. 

5. J. W. Craggs', "The Determination of Capacity for Two- Dimensional Systems of Cy- 

lindrical Conductors," Quart. J. Math. Oxford, Series 1,17 (1946), pp. 131-137. 

6. C. M. Miller, "Capacitances of a Shielded Balanced-Pair Transmission Line," B.S.T.J., 

51, No. 3 (March 1972), pp. 759-776. 

7. H. W. Friesen, unpublished work. 

8. W, T. Anderson, unpublished work. 

9. T. A. Lenahan, "Experimental Test of Propagation- Parameter Calculations for 

Shielded Balanced Pair Cables," B. S.T.J. , this issue, pp. 627-636. 

10. L. V. Ahlfors, Complex Analysis, New York: McGraw Hill, 1953, p. 147. 

11. G. Hsiao and R. C. MacCamy, "Solution of Boundary Value Problems by Integral 

Equations of the First Kind," SIAM Rev., 15, No. 4 (October 1973), pp. 687-705. 

12. K. Hoffman and R. Kunze, Linear Algebra, Englewood Cliffs, N.J.: Prentice-Hall, 

1961. 



UNIFORM CABLES— II 625 



